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We show that the main features of the cuprates superconductors phase diagram can be derived 
considering the disorder as a key property of these materials. Our basic point is that the high 
pseudogap line is an onset of phase separation which generates compounds made up of regions with 
distinct doping levels. We calculate how this continuous temperature dependent phase separation 
process occurs in high critical temperature superconductors (HTSC) using the Cahn-Hilliard ap- 
proach, originally applied to study alloys. Since the level of phase separation varies for different 
cuprates, it is possible that different systems with average doping level p m exhibit different de- 
grees of charge and spin segregation. Calculations on inhomogeneous charge distributions in form 
of stripes in finite clusters performed by the Bogoliubov-deGennes superconducting approach yield 
good agreement to the pseudogap temperature T*(p m ), the onset of local pairing amplitudes with 
phase locked and concomitantly, how they develop at low temperatures into the superconducting 
phase at T c (p m ) by percolation. 

PACS numbers: 74.72.-h, 74.80.-g, 74.20.De, 02.70.Bf 



I. INTRODUCTION 



After almost twenty years of research, the high criti- 
cal temperature superconductors (HTSC) still remains as 
an unsolved problem 0, 0, 0- All HTSC have a similar 
universal complex electronic phase diagram: the parent 
(undoped) compound is an antiferromagnetic Mott insu- 
lator, the superconducting phase has a dome shape at 
low doping and low temperature. The normal phase has 
a pseudogap in the underdoped region at low tempera- 
tures and a metallic phase at high temperatures in the 
overdoped region. Understanding the complexity of such 
normal phase is believed to be essential for solving the 
mechanism of HTSCQ. 

Another intriguing fact is the question of the inhomo- 
geneities in these materials; some families of compounds 
have a high inhomogeneous electronic structure which 
display either stripe[J], patchwork|g, checkerboard[6j, or 
other forms 0- On the other hand, depending on the type 
of experiment, there are some HTSC materials that ap- 
pear to be more homogeneous or, at least do not display 
any gross inhomogeneityd, 0. It is possible that these 
distinct features are due to a phase separation transi- 
tion that produces different degrees of local hole doping 
densities and consequently, different properties. We have 
already discussed this possibility in a previous paper [To) 
and we now perform more detailed calculations in con- 
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nections with several new experimental data. 

Recent angle resolved photoemission (ARPES) ex- 
periments with improved energy and momentum 
resolution^, 0, 0, have distinguished two elec- 
tronic components in fc-space associated with the 
La,2- x Sr x CuOi (LSCO) system: a metallic quasi par- 
ticle spectral weight at the (tt/2,tt/2) nodal direction 
which increases with hole doping and an insulator like 
spectral weight at the end of the Brillouin zone straight 
segments in the (tt, 0) and (0,7r) antinodal regions which 
are almost incentive to the doping level. Compari- 
son with the non superconducting La,2- x -yNd y Sr x Cu04 
system^^? m which static stripes were first observed^, 
demonstrated that the antinodal spectral weight behav- 
ior is compatible with a quasi one dimensional elec- 
tronic structure where the hole rich stripes behaves 
as one dimensional metals and the hole poor stripes 
as insulators. These features demonstrated that in 
these compounds there are two aspects of the electronic 
structure flU H^ . [Tip. IT^| . Moreover, the large shift of the 
ARPES spectra[llj at the Fermi energy, what is called 
the leading edge shift and it is interpreted as the super- 
conducting gap is maximum at the antinodal region. This 
is an indication of the d-wave symmetry of the supercon- 
ducting order parameter, the zero temperature supercon- 
ducting gap Ao vanishes in the nodal but it is maximum 
in the antinodal directions. 

Another technique which has been refined and revealed 
new aspects of HTSC is scanning tunneling microscopy 
(STM). It is complementary to ARPES because it probes 
the differential conductance or the pairing amplitude A 
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directly on the surface of the compound. New STM 
data with great resolution have also revealed strong in- 
homogeneities in the form of a patchwork of (nanoscale) 
local spatial variations in the density of states which 
is related to the local superconducting gap0 0, Ht| . 
More recently it was possible to distinguish two dis- 
tinct behavior: well defined coherent and ill-defined in- 
coherent peaks depending on the exactly spectra loca- 
tion at a Bi 2 Sr2CaCu 2 8 +s (Bi2212) surface H3 EE El- 
STM experiments have also detected a regular low en- 
ergy checkerboard order in the electronic structure of the 
Bi2212 family above the superconductin g c ritical tem- 
perature (Tc)|2l| and at low temperaturepof and in the 
Na x Ca 2 - x Cu0 2 Cl 2 ^. 

A third important set of experiments to describe the 
HTSC phase diagram is the tunneling current p3. |24| . 
New techniques have recently shown the existence of two 
energy gaps which behaved differently under an applied 
magnetic field |25l. l26l l27 | . Tunneling experiments using 
superconductor insulator superconductor (SIS) with in- 
sulator layers with various size and resistivity have also 
shown distinct sets of energy scales and have also led to 
the idea that the richness of the phase diagram as func- 
tion of doping may be due to charge inhomogeneities and 
charge cluster of different size in the Cu-0 planes [23, |2!j. 

The ARPES and STM experiments, are surface probes 
what may suggest that the inhomogeneities may be a 
surface effect, but charge disorders were also detected 
by bulk experiments, like the stripe phase in materi- 
als similar to LSCO|3, H3 and also local variations in 
the charge [3l| . Another bulk sensitive experiment is 
nuclear magnetic and quadrupole resonance (NMR and 
NQR) which have provided ample evidence for spatial 
charge inhomogeneity in the Cu0 2 planes [12 HE l34| . 
Similarly, Singer at al|3i| measured a distribution of T\ 
over the Cu NQR spectrum in bulk LSCO which can 
be attributed to a distribution of holes p with a half 
width of Ap/p w 0.5. More recently, NMR results on 
Lai,g- x Euo. 2 SrCu04 were interpreted as evidence for a 
spatially inhomogeneous charge distribution in a system 
where the spin fluctuations are suppressed |35j. This new 
result is a strong indication that the charge disorder may 
be due to a phase separation transition. 

These unusual features of cuprates led to theoreti- 
cal proposals that phase separation is essential to un- 
derstand their physics [3^. l37| . In fact, phase segrega- 
tion has been observed on the La 2 CuOi+s by x-ray and 
transport measurements piU l39| |. They have measured a 
spinodal phase segregation into an oxygen-rich (or hole- 
rich) metallic phase and an oxygen-poor antiferromag- 
netic phase above T=220K. Below this temperature the 
mobility of the interstitial oxygen becomes too low for 
a further segregation. La 2 Cu04 + s is the only system 
where ion diffusion has been firmly established, although 
there are evidence of ion diffusion at room temperature 
in micro crystals of the Bi2212 superconductors at a very 
slow rate|40|. 

On the other hand, recent NMR studies on 



Y Ba 2 Cu^OQ +y (YBCO) have demonstrated a complete 
absence of static phase separation or at least an absence 
of gross inhomogeneities [a 0- Contrary to what was 
measured in Bi2212 and in LSCO, the maximum hole 
doping variation Ap found in YBCO is very smallQ- 
Loram et alQ analyzed also the specific heat of YBCO 
and Bi2212 and concluded that there are evidences for 
an uniform doping density in these materials. Surface 
measurements are difficult to be performed in YBCO 
and therefore they are not as conclusive as in Bi2212 
or LSCOHEl. 

In this paper we develop the idea of two distinct [4l| 
pseudogaps in which the lower one is associated with the 
onset of superconductivity We then take the upper 
pseudogap line as a line of phase separation transition, 
and introduce a model to make quantitati ve p redictions 
through the Cahn-Hilliard (CH) approach[ia El in 
section II. In this way, gross and weak disordered systems 
differ only by the degree of mobility or diffusion of the 
particles and different systems can belong to the same 
universality class. In section III we use the Bogoliubov- 
deGennes (BdG) local method to the superconducting 
problem to calculate the local pairing amplitude in meso- 
scopic clusters with random, stripes, checkerboard, Gaus- 
sian and other forms of inhomogeneities in the charge 
density. The results on stripe-like formations applied to 
the LSCO system reveal a pseudogap phase characterized 
by the building up of superconducting islands or puddles 
with, as in BCS theory, the phase locked. Consequently 
the superconducting phase is reached at low tempera- 
tures by the percolation of these islands [iM 14(1 147^ . The 
details of these calculations are discussed in section II 
and III with their consequences to HTSC in section IV. 



II. THE CH APPROACH 

Our main assumption is that the high pseudogap line, 
which we call T ps (p m ), falls to zero near p ss 0.20, and 
is independent of the superconducting phased ESj * s 
the onset of the phase separation. Timusk and Sttat 
called this line a crossover boundary^] and it is distinct 
of the lower pseudogap Thus a given system start 
to phase separate at T ps (p m ) and this process increases 
continuously as the temperature goes down. 

Thus the phase separation problem in HTSC is a dy- 
namical and depends strongly on the initial conditions, 
on the temperature and how the system is quenched be- 
low the phase separation line, the mobility of hole and 
ions, and so on. However, the information on most of 
these procedures are not available and on has to work 
out the phase separation process backwards, to wit, use 
parameters which yield the final configurations of stripe, 
checkerboard, or other patterns. An appropriate frame- 
work to st udy such process mathematically is by the 
CH the orv|43j . which we have already applied to the 
cuprates jlOj. 

The CH approach to phase separation was conceived to 
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describe the continuous transition of binary alloys, but, 
in principle, can be applied to any system that undergoes 
this type of continuous transition ji^. As we can infer 
from the stripe phases, in a compound with an average 
of p m = 1/8 hole per copper atom, the antiferromagnetic 
insulating phase has stripes of nearly zero holes per cop- 
per atom and metallic ones with larger values of local 
charge density. As we discuss below, this behavior can 
be well described by the CH theory. 

Starting with small fluctuations of the local charge 
density, the CH non-linear differential equation which 
describes the phase separation at a temperature T be- 
low the phase separation transition at T ps can be written 

as[T3,li3 : 

^ = -MV 2 (e 2 V 2 u + A 2 (T)u-B 2 u 3 ), (1) 
ot 

where u is the order parameter associated with the local 
variation p(x) in the average number of holes per copper 
atom p m at a given point x, defined as u(x) = p(x) — p rn 
and we expected u(x) » above and near the T ps . e, 
and B are fixed parameter, the parameter A(T) depends 
on the temperature T and B and the ratio ±A(T)/B 
yields the two equilibrium densities. M is the mobil- 
ity of the particles and dictates the time scale of the 
phase separation process. Compounds with larger val- 
ues of M phase separate easier than those with smaller 
values, and, therefore, it will differentiate among the dif- 
ferent HTSC families. As the temperature goes down 
below T ps (p m ), the two equilibrium order parameter (or 
densities) spread apart from one another and the energy 
barrier between the two equilibrium phases E g (T,p m ) 
also increases. E g — A 4 (T)/B which is proportional to 
(T ps — T) 2 . E„ can be identified with the upper pseu- 
dogap signal [l0[ which is detected by several different 
experiments ja. l42fl . A discussion on the details of the con- 
vergence criteria, boundary conditions and dependence 
on the initial parameters in one, two or three dimensions, 
can be found in our previous work [44). 

In Figs.Q) an d @ we display the mapping of the or- 
der parameter for a system with p m = 1/8 and with a 
maximum Ap = 1/8. Both figures have a small initial 
density fluctuation of Ap = 0.02 but Fig.Q has more 
random initial conditions than the more symmetric case 
of Fig.©. We notice that, for larger times they are al- 
most indistinguishable but for short times they are very 
different. It is interesting that, despite the large time 
evolution, the system keeps very symmetric patterns. In 
Fig. (J2Jl) we can see a mixture of checkerboard and stripe 
formations which tends to evolve into a pure stripe phase 
(Fig. (2b)). As the phase separation process continues, 
the systems tends to a complete phase separation and 
the early symmetric patterns are lost. This degree of 
charge phase separation can be reached very fast if the 
holes have a large mobility M. This finding enable us to 
speculate that the differences measured among the HTSC 
families are due to different values of M. By the same to- 
ken, the differences found in compounds of a given family 



can be attributed to differences in the quenched process, 
that is, the rate of cooling below T ps (p). 




FIG. 1: (color online) The mapping of the density order 
parameter during the process of phase separation. The start 
order parameter is small random variations around u — 0. 
The order of the figures is a) t=400 time steps, b) t=800, c) 
t=1000 and d) t=4000. 




FIG. 2: (color online) The same time evolution of Fig.l but 
with a more symmetric initial condition around u = 0. The 
order of the figures is a) t=400 time steps, b) t=1000, c) 
t=4000 and d) t=20000. For shorter times, the phase sepa- 
ration process develops symmetry patterns which are lost for 
larger times. 

Another possible way to follow the whole phase separa- 
tion process is to analyze the charge histogram evolution 
with time. Thus, in Fig.ljJJ we show the phase separation 
progress in terms of the histograms of the order parame- 
ter. One can see the tendency, as the time flows, to form 
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a density pattern of a bimodal distribution around the 
two equilibria conditions p(i)± = ±A/B. Thus, systems 
with a high mobility will probably reach a state where 
there are two type of regions with high and low densi- 
ties. This is similar to the stripe phase in the LSCO 
system. 
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H = - ]T t io c\ a c ]a + ^(//i)n jCT 

+U n^nn + n icr n ja/ (2) 

i (ij)crcr' 

where c\ a (ci a ) is the usual fermionic creation (annihi- 
lation) operators at site Xj, spin cr{||}, and rii a = 
c ia .Ci a . tij is the hopping between site i and j. Here 
we have implemented in our calculations hopping val- 
ues up to 5 th neighbors derived from the ARPES data 
for YBCO 53]. In their notation, the hopping param- 
eters are: t = ti=0.225eV, i 2 /ti=-0.70, i 3 /ti=G.25, 
i 4 /*i=0.08, t 5 /ti=-0.08. U = l.lt is the on-site and 
V = — 0.6t is the nearest neighbor phenomcnological in- 
teractions. fXi is the local chemical potential. All the 
calculations presented here use the same set of parame- 
ters, and clusters with periodic boundary conditions. 



FIG. 3: The histogram time evolution of the density order 
parameter of Fig.l with a symmetric initial condition around 
u = 0. 

Therefore, taking the large pseudoeappil ITol Esj as the 
phase separation temperature T ps , it is possible to infer 
that the different HTSC unfold in different patterns as 
described by Figs. 1 and 2. As a consequence, since un- 
derdoped compounds have a very high T ps , they phase 
separate into a bimodal distribution faster than the op- 
timally and overdoped compounds. The NQR results of 
Singer et ah33] suggest a bimodal distribution of charge 
in the LSCO system. 

Thus, in order to calculate the local critical tempera- 
ture for such an inhomogeneous system, we need to use 
the Bogoliubov-deGennes theory. This will be outlined 
in the next section. 



III. THE LOCAL SUPERCONDUCTING 
CALCULATIONS 

Here we discuss the main points of a local supercon- 
ducting calculation to deal with the effect of the charge 
disorder which follows directly from the CH patterns de- 
scribed above. The general way to perform this, in a 
system without spatial invariance, is through the BdG 
mean-field theory which has been largely used in the 
HTSC problemHaESilil. The important and novel 
point introduced here is that we take the initial charge 
distribution derived from the CH results. The procedure 
starts with the extend Hubbard Hamiltonian 




U„(Xi)\ / U n (Xi)' 

\=En\ 

V n (Xi) I \.V n (xi) 



(3) 



These BdG equations are solved self-consistcntly to- 
gether with the pairing amplitude |49| 



Aa(xi) = -U ^2 u n (Xj)v* (Xj ) tanh 



2k B T' 



(4) 



y r ^ 

Aa(xi) = -— 2_Jit„(xi)v*(xi + 5) + v* n (xi)u„(xi 

n 



+,5)]tanh 2fc B T' 
and the hole density is given by 



(5) 



1 



2^[| W „(xO| 2 /„ + |«n(x. ( )| 2 (l -/„)], (6) 



where f n is the Fermi function. 

We have performed self consistent calculations with 
Eqs.JSJ) and JBJ on clusters up to 24 x 24 sites with ho- 
mogeneous and inhomogeneous local doping. The major 
difference from previous BdG calculations is that instead 
of an impurity potentially 0H1 to account for the charge 
disorder, we have fixed the initial local charge densities 
throughout the calculations in order to take into account 
the results of the CH approach and also in agreement 
with the data on stripe formation0,[3(]]. This novel pro- 
cedure is necessary to study the formation of the super- 
conducting regions on the local density patterns which 
results from the phase separation process as shown, for 
instance, in Figs.l and 2. In the next section, we perform 
the superconducting calculations with d-wave symmetry 
for clusters with disordered local charge. 
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IV. RESULTS 

We start with BdG calculations on clusters of uniform 
density ranging from zero to p m = 0.3 with the hopping 
parameters which make the values of T c (p m ) vanish near 
these limits as listed below Eq.2. Then we perform cal- 
culations with inhomogeneous clusters and concentrate 
on the stripe geometry which occur in LSCO0. We con- 
nect the stripe phase with the above CH results through 
a scheme displayed in Fig.© where the value of X is re- 
lated with the degree of the phase separation. Due to the 
large values of T ps {p m ) at the strong underdoped regime, 
the level of phase separation is maximum at low temper- 
atures, and X — for doping values of p m < 0.05. For 
0-12 < p m < 0.19 the value of X may build up to 0.05. 
Notice that when X=0, the charge distribution is a bi- 
modal and the system is divided into two distinct regions; 
insulator and metallic. In this case the metallic regions 
are in the limit of percolation but for X=0.04 or big- 
ger, despite of the insulator regions present, the metallic 
character dominates over the entire system. This scheme 
is an approximate way to deal with the phase separation 
which leads to the stripe charge configuration of real sys- 
tems. 
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FIG. 4: The low temperature stripe profile of a 14 x 14 cluster 
with 196 sites used to model a compound of average doping 
level of p m . The BdG superconducting calculations are made 
on these clusters and < X < 0.06 measure deviations from 
a bimodal charge distribution. 



Thus, the goal is to study the local pairing gap at a 
site Xi or simply "i" as function of temperature A(i,T) 
on clusters with charge stripes. Following the values 
of A(i,T) it is possible to draw many interesting con- 
sequences to the phase diagram. To explain our ap- 
proach, we will firstly analyze the results for a cluster 
with p m = 0.16inal4xl4 cluster as shown in Fig. © ■ At 
high temperatures but below T ps (p m = 0.16) « 400K[3|, 
the system is a disordered metal going through a continu- 
ous phase separation as the temperature is decreased. At 
low temperatures (< 200-ftT) this compound may be com- 
posed of six stripes with local charge density p(i) given 
by 0.0 - 0.05 - 0.0 - 0.16 - (0.32 - 0.05) - 0.16. The low 
density regions are at the left and the high density are 
at the right of the cluster and one follows the other by 
periodic boundary conditions. These two markedly dif- 
ferent regions were clearly detected by several ARPES 
measurements 0, 0, 0, 0|. At high temperatures 
there is no superconducting regions in the sample, but 
at T* = 65K, as we can see in Fig.©, some local super- 
conducting gap arises forming superconducting islands at 
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FIG. 5: (color online) The temperature evolution of the local 
pairing gap amplitude A(i,T) (in units of eV) for systems 
with stripe disorder and average doping level of p m — 0.16 
with T* « 65K and T c « 42K. 



the most dense or metallic regions. Thus this tempera- 
ture is the onset of pair formation or the beginning of 
the pseudogap phase but these superconducting regions 
are isolated from one another in the disordered metallic 
matrix lit! l47j . Upon cooling down, the supercon- 
ducting regions become more robust and new ones are 
built up at the metallic and also where there is a hole 
fluctuation (at X in Fig.©) at the lightly doped regions. 
The induced pairing amplitude at the insulator regions 
is an unexpected phenomena and it is remarkable that 
it is the origin of the superconducting phase through the 
percolation of the local superconducting order parameter 
at different locations with their phase locked. As seem 
in Fig.©, around 42K the pairing amplitudes develop 
also at X = p(i) = 0.05 but for the p m = 0.06 sam- 
ple, they develop even at X = p(i) = 0.001 (see Fig.©) 
and thus, the superconducting regions cover more than 
50% of the CuOi plane. Consequently they percolate [54^ 
and the whole system is able to hold a current without 
dissipation, that is, T c (p m = 0.16) = 42K is the super- 
conducting critical temperature for this compound. 

Notice that the assumption that the pairing ampli- 
tudes have a rigid phase, as in a BCS superconductor, 
has also experimental support [5^, although it is against 
the phase disordered scenario which lately has gained in- 
creased theoretical interest fH IHril] . 

We now apply this analysis to a series of compounds in 
order to show how the main features of the whole LSCO 
phase diagram can be derived. Our results are shown in 
Fig.© for mostly underdoped samples and in Fig.© for 
larger average doping values. For very lightly doped com- 
pounds like p m = 0.03, due to the high values of T ps (p m ), 
there are stripes of only p{i) — 0.0 and p{i) — 0.06 sep- 
arated by a small boundary with p(i) = 0.03. We label 
a region as metallic if it has a density of p(i) > 0.05, 
since the doped regions with p(i) = 0.06 occupy less 
than half system, what is lower than the two dimen- 
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sional percolation limit of 50-60%[54j, this compound 
is not a metal, although it has a metallic behavior at 
high temperatures Such behavior can be explained 
as due to the holes which can tunnel over the dense 
p(i) = 0.06 regions. This tunneling can also be the ori- 
gin of the zero temperature pseudo gap (ZTPG) detected 
by STM plll9ll20| . For these underdoped compounds, 
the pairing amplitude develops strictly in these metallic 
or heavy doped regions (see top panel of Fig.JTJl), and 
the superconducting islands occupy less than 50% of the 
available area, what is below the percolation limit. 
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FIG. 6: (color online) The temperature evolution of the local 
pairing amplitude A(i, T) (in units of eV) for systems with 
stripe charges with average doping level ranging from p m — 
0.03 to p m = 0.12. The onset temperature (T c ) of percolation 
is shown in every panel 



For a compound with p„ 



0.06, T ps is still very 



high, the density profile is characterized by a very small 
charge fluctuations around a bimodal distribution given 
by X = 0.001. The stripe regions have densities given 
by p(i) = 0.0 - 0.001 - 0.0 - 0.06 - 0.119 - 0.06. This 
system can be considered in the metallic limit which is 
50% of the sites with p(i) > 0.04. We see that the on- 
set of superconducting islands is at T* = 100K and the 
percolation threshold occurs at T c — 15K, specially due 
to the already mentioned unexpected pairing amplitudes 
induced at the low density sites with p(i) — 0.001 (see 
Fig.©). Thus for temperatures between T — 15K and 
T = 100K the system is a poor metal with insulator, 
metallic and superconducting regions. The presence of 
these superconducting islands in many compounds are 
verified by several different experiments. Perhaps the 
most clear demonstration of these static superconduct- 
ing regions is through the tunneling data |22 , l29l . More 
recently, measurements of the Nernst effect 58j demon- 
strated also the presence of the local superconducting re- 
gions above T c {p m ) although it was interpreted as due to 
superconducting fluctuations instead of the static cluster 



considered here. For compounds with an average doping 
larger than p m = 0.06, we notice that the onset of su- 
perconductivity T* decreases almost steadily, while the 
onset of percolation T c goes through a maximum at the 
optimum doping p m = 0.16. The reason for this behav- 
ior in our calculation is the deviations from the bimodal 
distribution given by the increase of the phase separation 
parameter X. 
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FIG. 7: (color online) The temperature evolution of the local 
pairing gap amplitude A(i,T) (in units of eV) for systems 
with stripe disorder with average doping level ranging from 
p m = 0.14 to p m = 0.22. Notice that the p m = 0.22 com- 
pound does not undergoes a phase separation transition and 
its density fluctuates around p m = 0.22. 

In the Fig.Q we can see that the phase separation 
process and the local superconducting calculations fol- 
low a similar pattern up to p m — 0.20. Following current 
trends |M l42l |48| . the phase separation ends at p m « 0.20 
and, for heavier doped compounds, the charge disorder 
is very weak like a small fluctuation around the aver- 
age value p m , similar to the charge distribution shown 
in Fig. (3a). This is the same kind of charge fluctuation 
which occurs above the phase separation transition tem- 
perature T ps . Consequently, the compound with p m — 
0.22 differs greatly from the others compounds shown in 
Fig.Q, and it has only local densities p(i) w 0.22 which 
is in the metallic range. The pairing gaps A(i, T) are 
built in the whole system at T — 30K and not in islands 
or droplets. If we adhere to the assumption that the 
T* is the onset of superconducting correlations, we see 
that, for these p m > 0.20 overdoped samples T* merge 
into T c . The other important consequence is that the 
normal phase is much more homogeneous, without insu- 
lating and superconducting regions, than the p m < 0.20 
compounds. This is seem experimentally through the 
Fermi liquid behavior of many measurements carried in 
the heavy overdoped regimej^Q. 

In Fig.© we used the calculated local A(i,T) on each 
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p (average doping) p(doping level) 



FIG. 8: (color online) The onset of superconducting islands 
temperature T*(p m ) and the percolation temperature T c (p m ) 
taken from Figs.JJJ and @. It is also shown the phase sepa- 
ration (dashed) line or upper pseudogap pi. EUl . 



compound to derive the LSCO phase diagram, to wit, the 
onset of superconducting temperature T*{p m ) (squares) 
and the percolation temperature T c {p m ) (circles) as func- 
tion of p m , as derived from Figs.Q) and ©. The values of 
T*(p m ) are in good agreement with the measurements at- 
tributed to the lower pseudogap which is usually related 
with the onset of a superconducting p roperty 13, as, 
for instance, the tunneling results |22ll29| and the Nernst 
effect |58|. The values of T c (p m ) is also in good quanti- 
tative agreement with the experimental superconducting 
phase boundary pi l3l l58| . 

We have also taken the larger superconducting gaps at 
low temperatures for each compound from Figs. I© and 
(Q) and plotted in Fig. © . The calculated maximum gap 
for each compound at low temperature (Arj(p m )) is in 
reasonable agreement with the ARPES zero temperature 
leading edg e shift or the magnitude of the superconduct- 
ing gap 14] . In this figure, just to show the effect of the 
disorder in our calculations, we have shown the values of 
Ao(Pm) for homogeneous compounds. We see that the 
disorder increases the average zero temperature gap dra- 
matically at low doping and in a weaker way, in the far 
overdoped regions. The discrepancies around optimally 
doped samples between the experimental values and our 
calculations are likely to be due to our approximate stripe 
configurations on small clusters. 

It is important to notice that this stripe scheme is able 
to capture the very curious dual behavior of the electronic 
structure in LSCO systems [12j. At the lightly doped 
regime, due to the high values of the phase separation 
energy barrier E g (T,p m ), the charges move preferably 
along the high densities stripes, and the Cooper pairs 
are formed along them, as demonstrated in Figs.© and 
O by the calculated A(i,T — OK). Consequently, with 



FIG. 9: The zero temperature superconductor gap as func- 
tion of the doping level. The diamond points are for a cluster 
with uniform density. The square are for an inhomogeneous 
cluster according the CH results. The circles are the experi- 
mental leading edge gap from Ref.(IT^|'l. 



a k-space probe, the superconducting gaps for lightly 
doped samples are measured mainly in the (tt, 0) and 
(0, 7r) antinodal regions and the spectral weight segments 
are entirely near these antinodal regions. This behavior 
is expected for ID stripes 12] and the measured values of 
the zero temperature leading edge as function of p m can 
be qualitatively reproduced by the An(p m ), as shown in 
Fig.©. 

Since by assumption E g (T,p m ) is closely related with 
the high pseudogap and, therefore decreases rapidly with 
p m , it is easier to the holes in compounds near optimally 
doped than underdoped ones to tunnel over among the 
high and low doping stripes, yielding a 2D character to 
these systems. Thus as p m increases, the samples tend to 
change continuously from ID to a 2D metallic behavior 
and this is detected by the increase of the spectral weight 
near the Fermi surface along the [1,1] nodal direction |12|. 
Another evidence of the charge tunneling between dif- 
ferent stripes is the presence of the incoherent ZTPG 
measured on the surface of Bi2212 compounds [l8l ITs| 
which should scaled with the energy barrier E g (p m , T) . 
The ZTPG are more frequently than the superconduct- 
ing gap in the underdoped region and disappear near 
p m — 0.190, |23. Evidence of these two types of 
gaps in HTSC compounds came also from tunneling ex- 
periments carried with different resistances [281 129| . 



V. CONCLUSIONS 

We have worked out a complete scenario to HTSC and 
provided an interpretation to the upper and lower pseu- 



8 



dogap line and also to the superconducting phase. Taking 
the upper pseudogap as the phase separation temper- 
ature, we have calculated the local pairing amplitudes 
A(i, T) which, as in BCS, is assumed to have their phase 
locked and the superconducting phase is reached by per- 
colation at T c (p m ). In this way, we derived the phase di- 
agram of the onset of A(i, T), that is T*(p m ) and T c (p m ) 
for a HTSC family. Although the process of phase sepa- 
ration varies continuously with the temperature and de- 
pends on the sample preparation, for simplicity, we use in 
our calculations only the low temperature static config- 
uration. Despite this simplification, the method is quite 
general, and to demonstrate it, we reproduced results in 
good agreement with the LSCO series. 

The many values of A(i,T) in a single compound 
agrees with the several recent very fine STM maps on 
Bi2212 surfaces. The calculations with the CH stripe 
configuration have also provided novel interpretation to 
important aspects involved in the electronic structure 
of this type of disorder: the inhomogeneous electronic 
dual nature of underdoped LSCO cuprates measured by 
the straight (ID) segments near the anti nodal regions 
and the spectral weight near the nodal regions where the 
Fermi surface develops. The presence of the ZTPG peaks 



in several Bi2212 compounds, as measured by high reso- 
lution STM data, may be due to the energy barrier be- 
tween the low and high density regions in a given sample 
and therefore, in this view, is connected to the phase sep- 
aration and not with the superconductivity. Our results 
indicate that the normal phase of cuprates is a disordered 
metal for p m < 0.2 composed by the coexistence of in- 
sulator and metallic regions below the phase separation 
temperature T ps . These regions are composed of non 
constant hole density regions which yields a variety of 
local T c (p(i)) at low temperatures. This inhomogeneity 
is the source of scattering between these non-uniform re- 
gions, Andreev reflection, etc, what is the cause of many 
nonconventional transport properties@,0- The study of 
these normal phase properties will be matter of future 
publication. 

Thus, in short, with the combination of the CH phase 
separation model with the BdG superconducting theory, 
we have shown that the phase diagram and some non- 
conventional properties of HTSC receive a coherent inter- 
pretation based on a temperature dependent static phase 
separation. 

This work has been partially supported by CAPES, 
CNPq and CNPq-Faperj Pronex E-26/171. 168/2003. 
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